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Abstract. We classify all stationary axi-symmetric solutions of topologically massive 
gravity into Einstein, Schrodinger, warped and generic solutions. We construct explicitly 
all local solutions in the first three sectors and present an algorithm for the numerical 
construction of all local solutions in the generic sector. The only input for this algorithm 
is the value of one constant of motion if the solution has an analytic centre, and three 
constants of motion otherwise. We present several examples, including soliton solutions 
that asymptote to warped AdS. 
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1. Introduction 

Topologically massive gravity |TJ|2) (TMG) is a 3-dimensional theory of gravity that exhibits 
gravitons and, for negative cosmological constant, black holes |3j. Its action is given by 

s TMG = -J d3x V\9\{ R +^)-^ 1 j ^ xe ^ u T ^p (V p - + \ r V rT -) (!) 

The interest in TMG was rekindled in 2008 due to a paper by Li, Song and Strominger j4], 
and the discussion that ensued after the papers by Carlip, Deser, Waldron and Wise [5] and 
two of the current authors j6]. 

In this work we continue the TMG saga by explaining how to construct all its stationary, 
axi-symmetric solutions. Finding non-trivial solutions to the TMG equations of motion is 
a very difficult problem and thus far few solutions are known. For a recent summary 
and algebraic classification of all known exact solutions see [F| and references therein. As 
remarked in one important reason for the scarceness of known solutions is the presence of 
no-go theorems. For instance, there are no non-trivial static solutions. In fact all solutions 
that have a (non-null) hypersurface-orthogonal Killing vector are Einstein (§]. 

A weaker assumption that can aid in finding solutions is to require instead two 
commuting Killing vectors. Such spacetimes are the focus of this paper. In this case a 
powerful method pioneered by Clement |9] becomes available: stationary axi-symmetric 
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TMG can be formulated in terms of a 0+1 dimensional theory. We call this theory 
"topologically massive mechanics" (TMM). Despite its superficial simplicity it is surprisingly 
difficult to find analytic solutions of this theory. On the other hand there are some non- 
existence results: for the tuning \fi£\ = 1 all solutions that satisfy Brown-Henneaux boundary 
conditions are Einstein [To] and for the case \fj£\ = 3 all solutions that satisfy a set of 
asymptotically warped boundary conditions are the null-warped black holes (there is an 
anlog statement for fx£ > 3 asymptotic Schrodinger black holes) [IT]. Let us classify the 
solutions of TMM into four sectors, restricting for the time being to negative cosmological 
constant i 2 > 0. 

(i) Einstein: these solutions solve also the 3-dimensional Einstein equations; the only 
solutions in this sector are locally AdS 3 and obey Brown-Henneaux boundary conditions 

(ii) Schrodinger: these solutions are either asymptotically Schrodinger or asymptotically 
AdS 3 . They behave like zero modes of pp- waves. 

(iii) Warped: these solutions are locally or asymptotically warped AdS3 

(iv) Generic: all solutions that are neither Einstein, nor Schrodinger, nor warped 

The three first categories were essentially exhausted in [9] and correspond to all presently 
known solutions of TMM. 

In this work we analyse TMM and its space of solutions by several methods. First we 
perform a Hamiltonian analysis and count the dimension of the physical phase space. We 
find that the phase space in general is six-dimensional. However, we find an interesting 
four-dimensional sub-space of the physical phase space occurring naturally. This subspace 
is generated by two conditions that are preserved by the equations of motion. Assuming 
that these conditions hold reduces the dimension of the physical phase space to four. For 
this special subspace the Hamiltonian dynamics simplify considerably and we are able to 
construct all local solutions in closed form. In fact, we prove that one is left with exactly 
the three first sectors: Einstein, Schrodinger and warped. The generic sector then contains 
all solutions outside this sub-space. 

To describe the generic sector we resort to a numerical analysis. It turns out that one 
needs to employ three slightly different methods depending on the initial data. We provide a 
complete algorithm for numerically constructing a solution given any consistent initial data. 
The numerical solutions that we find in the general sector have non-constant curvature 
invariants. Some, but not all, of the examples we construct develop naked curvature 
singularities. In the set of smooth examples we discuss in particular soliton solutions that 
asymptote to warped AdS and unravel their asymptotic (non-analytic) behaviour. We also 
generalise many of our results to TMG with vanishing or positive cosmological constant and 
address possible applications to recent new massive gravity theories. 

This work is organised as follows: In Section [2] we review the derivation of TMM 
and provide a useful reformulation of the equations of motion. In Section [3] we perform 
a Hamiltonian analysis of TMM and count the dimension of the physical phase space. In 
Section [| we classify the solutions into Einstein, Schrodinger, warped and generic and prove 
that the generic sector is identical to the solutions violating the two conditions found in 
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the Hamiltonian analysis. We provide simple algorithms to construct local solutions for all 
possible cases in Section [5} In Section [6] we give several examples and show in particular 
that the generic sector is not empty. An outstanding class of examples are soliton solutions 
that asymptote to warped AdS. In Section [7] we provide generalisations to the case with 
vanishing or positive cosmological constant and give directions for future work. 

Before starting we list our conventions. We use signature (+,—,—). For the target 
space of TMM we bring the 3-dimensional Minkowski metric into light-cone gauge 2r] + _ = 
— fjYY = 1- We fix the sign of the target space epsilon-tensor by choosing e+_y = +|. 
To reduce clutter we drop overall factors in front of actions, since we are only interested in 
equations of motion and their solutions. Our sign convention for the Ricci-scalar is such that 
R < for AdS. When we write Lorentz vectors as row or column vectors the components of 
these vectors always have upper indices. We define the cross product between two Lorentz 
vector as (A x B)* = e l jkA j B k . 



2. Topologically massive mechanics 



Following Clement [9] (see also 10 12 13]) we employ the following gauge for a generic 
stationary, axi-symmetric 3-dimensional line-element 

e 2 

ds 2 = dx* dx u + dp 2 (2a) 

with the 2-dimensional metric 

9r=\ ; : ) 




whose determinant is given by the norm of the Lorentzian 3- vector X = (X + , X , Y) 

det g = X 2 = XXi = X i X j r} i j = X + X~ - Y 2 (2c) 

All quantities depend on the radial coordinate p only. The vector X must be spacelike 
everywhere, possibly except at boundaries of the range of definition of p where it can become 
null. These boundaries can correspond to curvature singularities, Killing horizons, the radial 
centre or the asymptotic boundary of spacetime. For brevity, we use the term "centre" to 
refer to any point where X 2 = 0. 

Plugging the Ansatz for the line-element ^ into the action ([T| and integrating out the 
Killing coordinates x M yields a reduced action in 0+1 dimensions: 



'TMM 



dp e ( \ e- 2 X 2 - | - i- e" 3 e ijk XPx k ) (3) 



Dots denote derivatives with respect to the radial coordinate p, which acts as "time" in 
the particle mechanics system described by the TMM action (|3]). Any solution to the field 
equations descending from the action ^ can be oxidised to a solution of TMG by virtue of 



All stationary axi- symmetric local solutions of topologically massive gravity 



4 



the line-element ([2]). All stationary axi-symmetric solutions of TMG can be obtained locally 
in this way (9). 

For [i — > oo the action ^ simplifies to the geodesic action. For finite \i this geodesic 
action is corrected by a non-linear higher derivative interaction originating from the Chern- 
Simons term. TMM describes a spacelike particle moving in an auxiliary 3-dimensional 
Minkowski target space. The Einbein e in the TMM action (|3]) acts as a Lagrange multiplier. 
We set it to unity after variation, without loss of generality. This amounts to a gauge fixing 
of the p-reparameterisation invariance. Variation with respect to the Einbein yields the 
Hamilton constraint 

G=-X 2 + --- e ijk X l Px k = (4) 
2 i l /i 

Variation with respect to the target space coordinates X % establishes the equations of motion 

*i = -j- £m (3M fc + 2 A^X fc ) (5) 

Invariance of the TMM action ^ under Lorentz transformations leads to the conserved 
angular momentum «/, which provides first integrals to the equations of motion 

J = XxX + -Xx(XxX) -ix(Ixi) (6) 

We provide now a useful reformulation of the equations of motion. The conserved angular 
momentum can be represented as 

J% = e ijk X?X k - i- (5X 2 + ^) Xi + ±- (XX) Xi-- X 2 X, (7) 

Combining the Hamilton constraint Q conveniently with the first integrals ^ yields the 
condition 

e ijk JX>X k = \ X 2 X 2 - (II) 2 - j 2 X 2 (8) 

Contracting the angular momentum J with X and using the equations of motion leads to 
the constraint 

XJ = — ({XX) 2 -X 2 X 2 ) (9) 

Given some initial data X, X and values for the constants of motion J compatible with 
the conditions ^ and (|9| the first integrals ^ provide an evolution equation for X. Note 
that no third [second] derivative terms are present in (|7|) [(J8|, (J9|]. 

The equation for the angular momentum ^ together with the constraints ([8]), (|9| fully 
describe the dynamics of the theory defined by the TMM action ^ at all points, except 
at singularities of the line-element ([2]). Our goal is to find all solutions of these equations, 
and thus locally all stationary, axi-symmetric solutions to TMG. We start with a counting 
of the physical degrees of freedom. 
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3. Hamiltonian analysis 



To determine the dimension of the physical phase space we perform a constraint analysis in 
a suitable Hamiltonian formulation. We rewrite the TMM action ^ in first order form: 



TMM 



H{e,x\y\z\pf, 



Pi) 



(10) 



Our phase space is 20-dimensional, with canonical coordinates (e, x l , y l , z l ) and canonical 
momenta (p e , p* , p\ , pf) . The quantity x l is identical to X\ while y l and z % for e = 1 
correspond to the velocity X and acceleration X, respectively. Similar to the canonical 
analysis of TMG at the logarithmic pointj^] \pl\ = 1 14 we have identified already the 



appropriate canonical momenta in order to avoid a larger phase space with more constraints, 



since such an approach eventually would reduce to the formulation (10) that we are using 
as starting point (see also [15) ). The momenta p* and p\ are in fact Lagrange multipliers 
enforcing the relations x l = ey l and y l = ez l . The canonical Hamiltonian is given by 



H = eG 



with the Hamilton constraint 



1 2 1 

G = p*y l + pf z l - - y iV l + ^ + 7^ tijk x l y J z k 



Defining the quantities 
1 



p y i + 



2/i 



^ijk ^ y 



Pi - Vi tijk X 3 Z k 



Zi + — e ijk y 3 z k 



we obtain the following set of constraints (~ means weakly vanishing): 

Primary: pi w Hamilton: G ~ 

Ternary: $j w Quaternary: ~ 



Primary: p e w 
Secondary: Hj ps 



(12a) 



(12b) 



(12c) 
(12d) 



Two linear combinations of these twelve constraints are first class: p e and the extended 
Hamiltonian H ext 

H cxt = e (G + \y t ) (13) 

where the Lagrange multipliers A* are determined functions of the canonical variables. All 
other constraints are second class. We collect in |Appendix A the non-vanishing Poisson 



brackets between the constraints (A.l), results for the Lagrange multipliers A* (A. 2) and 



the Hamilton equations of motion (A. 3 ). We conclude that the constraints eliminate fourteen 



degrees of freedom and we are left with a six-dimensional physical phase space. 

The physical phase space has an interesting subspace. Namely, suppose that for some 



configuration all constraints (|T2j) are fulfilled and in addition the following constraints hold 

z% = z l Zi w (14) 



at some value of the "time" p 

Assumptions: y l ^i = y l Zi ps 
1 This is also known as "chiral point" or "critical point" . 
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Then these fourteen constraints hold for all times (see Appendix A). If x l , y l , z l are linearly 



independent then the extended Hamiltonian (13) reduces to the canonical one (11). The 



subspace of the physical phase space defined by the additional constraints ( 14 ) is only four 



dimensional. Interestingly, all currently known stationary, axi-symmetric solutions to TMG 
belong to this subspace. 

We note that the angular momentum ^ takes a rather simple form in terms of canonical 
variables. 

•/, ~ (.'V' • U J l> llh ) (15) 



The angular momentum commutes with the canonical Hamiltonian (11), with the extended 
Hamiltonian (13), with the primary constraints p e ,p* and with the secondary constraints 
Hi. The constraints (12) and the Hamilton equations of motion (A.3) taken together are 



equivalent to the Lagrangian equations of motion ((7j)-(j9]). 

In the next section we shall freely switch between the Hamiltonian and the Lagrangian 
formulation, depending on which description is most adequate to illuminate certain points. 
To facilitate this switching we always employ the gauge e = 1. 



In this gauge X % 



x\ X 1 



Gauge fixing condition: e 

y l and X % — z % . 



(16) 



4. Classification of all stationary axi-symmetric solutions 

We consider here all solutions to the equations of motion (|4])-([5]) [or, equivalently, of (|7j)-([9j]. 
We classify the solutions into four sectors: Einstein, Schrodinger, warped and generic. 

To simplify certain expressions we shall often exploit Lorentz transformations, constant 
rescalings 

p^VLp X^tt^X X^X X^ttX (17) 

and constant shifts of p, all of which leave invariant the line-element ^ up to 
diffeomorphisms (and possibly up to changes of ranges of definition of spacetime 
coordinates). The rescalings (17) together with the SL(2, R) transformations corresponding 



to Lorentz boosts and rotations form the group GL(2, \ 



4-1. Einstein sector 

This sector encompasses all solutions with vanishing acceleration, X = 0. 

The geodesic equations of motion are supplemented by a simple constraint on the norm 
of the velocity 



X = 



X 



4 

'J? 



;is) 



The solution to the equations of motion (18) is given by 



X = X( ) p + X(2) 



(19) 
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with some constant vectors X( ),X( 2 ), where X 2 ^ = — Jy. These solutions constitute all 



stationary, axi-symmetric solutions of Einstein gravity, as pointed out first by Clement [16] . 
By exploiting Lorentz transformations, rescalings of p and shifts of p we can bring X into 
one of the following forms (a 7^ 0) 

X=(a,^,±jp+i) X = (a,0,±- £ p) (20) 
The corresponding conserved angular momentum is given by 

J = ^ (1 T Id, 1(1 ± 0) J=l(lT /tf, 0, 0) (21) 



At p = the coordinate vector (20) becomes light-like. Note that by boosting in the 
(X + X~)-plane it is possible to eliminate also the single remaining constant of motion a. 
Thus all solutions are locally equivalent. We keep a arbitrary since it makes it simpler to 



connect (20) with the BTZ metric. The angular momentum (21) becomes light-like at the 
logarithmic point \p,£\ = 1 for the first case and is always light-like for the second case; we 
shall see below that the latter corresponds to extremal BTZ black holes. 



We discuss now briefly how to obtain the BTZ metric. Inserting the left solution (20) 
into the line-element ^ one obtains 

ds 2 = a (dx + ) 2 + - (dx-) 2 ± (e 2r - + e- 2r a) dx + dx~ - £ 2 dr 2 (22) 
a a 

where p = £(e 2r /a + e~ 2r a =F 2)/4 with a = AG(LL) 1 ^ 2 /£. The arbitrary constant a 
parameterises the angular coordinate <fi = (x + /a£ 2 ± x~)/2, which has periodicity of 2n. 
Defining new coordinates by x~ = =pi>, x + = £ 2 a u and parameterising the constant of 
motion as 1/a = 4G£L yields 

ds 2 = -4G£(L du 2 + L dv 2 ) - (£ 2 e 2r + 16G 2 Lle~ 2r ) du dv - £ 2 dr 2 (23) 



This is one of the standard ways to represent the BTZ black hole metric 17 , where L, L 



are related to outer/inner horizon r± and mass/angular momentum m,j as follows: 

L= (r+ + r ~ )2 L = (r+ ~ r - ? m = L + L j=L-L (24) 

16G£ 16G£ J v ; 



The extremal case L = arises as a singular limit of the first solution (20) and is described 
adequately by the second solution pOj ). Global (Poincare patch) AdS arises in the limit 
j — > 0, m — > —1/SG (j — > 0, m — > 0). The results above are obviously in accordance with 



the local triviality of 3-dimensional Einstein gravity 18 . All solutions in the Einstein sector 
are locally and asymptotically AdS. 

We shall encounter solutions in the other sectors that are also asymptotically AdS. In 
TMG there are several possible behaviours for the subleading terms. A general Fefferman- 
Graham like expansion for the coordinate vector in the limit of large p takes the form 



X = X m p + o(p) 



(25) 
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In TMG any vector X compatible with the expansion (25) leads to a 3-dimensional 
spacetime that is asymptotically AdS. If the subleading terms o(p) are constant then 
the more restrictive Brown-Henneaux boundary conditions are fulfilled. (This eliminates 



certain asymptotically AdS solutions of TMG.) These boundary conditions hold for (19), in 



accordance with 19 . If o(p) is logarithmic in p then the logarithmic boundary conditions 



of 20 21 are required. We shall encounter such a solution below. If o(p) behaves like 



p 1 £ , with e > 0, then yet other boundary conditions apply 21 . We give examples of such 
behaviour when describing the Schrodinger sector below. 

4-2. Schrodinger sector 

This sector encompasses all solutions where the coordinate vector X, its velocity X and its 
acceleration X are linearly dependent and acceleration is non-vanishing, X ^ 0. 

If X and X are parallel then the equations of motion have no solution with I/O, 
so we can disregard this case. Thus, in order to get linear dependence between coordinate 
vector, velocity and acceleration we assume with no loss of generality 

X = A(p) X + B(p) X (26) 

Then all triple products between any triple of coordinate vectors or its time derivatives 
vanish, and we obtain from the Hamiltonian constraint and its time derivative the conditions 

X' = -j 2 XX = (27) 

Contracting the equations of motion ([5]) with the coordinate vector and its acceleration 
establishes 

XX = X 2 = (28) 



The second equations in (J27I) , (28) imply that the assumptions (14) are fulfilled. Contracting 



the angular momentum ^ with the acceleration yields 

XJ = (29) 

Since the angular momentum is constant the last condition implies that the acceleration 
must be proportional to some overall function times a constant vector that has vanishing 
inner product with the angular momentum. We exploit Lorentz transformations to bring 
the light-like acceleration into a simple form. 

X = C(p) (1, 0, 0) (30) 

If C vanishes we are back to the Einstein sector discussed above. If C is non- vanishing we 



can integrate the acceleration (30) and impose the Hamiltonian constraint to determine X: 



X = (C(p), 0, ±h (31) 
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Now we integrate the velocity (31), impose the remaining constraints derived above and 



exploit constant shifts in p so that the norm of X is non-vanishing for positive p: 

X=(C(p), 0, ±jp) (32) 

Plugging this vector into the equations of motion ^ yields a third order ordinary differential 
equation for C(p) that can be solved straightforwardly. Exploiting rescalings of p we obtain 
the final result (\pJH\ ^ 1) 

X = [sp^/ 2 + ap + b, 0, ±- p) (33) 

which depends on two integration constants a, b and a sign s = ±. Again, a can be fixed to 
an arbitrary number by further Lorentz transformations and b can be rescaled by a positive 



number. We keep them arbitrary for convenience. The result (33) leads to 3-dimensional 
spacetimes with asymptotic Schrodinger behaviour]^] 

, 9 , „, ( ±2 dx + dx~ — dr 2 n (dx + ) 2 \ . . 

d * 2 Lo - * ( -2 + —IT-) ( 34 ) 



where /3 is a constant whose explicit value will be given below. The scaling exponent z is 
related to the parameters in the action as follows. 

* = ^ (35) 

For p£ = =1=3 we obtain the special case of (null-)warped AdS, z — 2. Despite the name 
"null warped" this solution belongs to the Schrodinger sector according to our classification 
since X, X and X are linearly dependent. The constraint ^ implies a light-like angular 
momentum 

J=(^(lTld)> 0, 0) (36) 

The angular momentum provides one constant of motion b. The periodicity of the angular 
coordinate <fi constitutes a second free parameter. Thus, again all geometries in this sector 
are parameterised by two constants, mass and angular momentum. 

Let us now obtain the Schrodinger line-element in standard form. We take the solution 
for X (33), plug it into the definition of the line-element ^ and obtain 



ds 2 = (sp^l 2 + ap + b) (dx + ) 2 ± ^ dx + dx~ - (37) 

The linear coordinate transformation x + = a(t ± £<f)), x~ = (±t — £(/))/ (2£a) together with 
the choice a = and the definition m T : = -2a 2 {sp^^' 2 + b) leads to 



'% - ^) dt 2 =f £m T dt d0 - (2p + i 2 ^) d0 2 - ^ (38) 
The Schrodinger solutions were discussed already by Clement [9] , except for the two logarithmic solutions 



( 39 1 , ( 40 1 below that were discovered in |22||23| . Higher-dimensional versions of the asymptotic Schrodinger 



line-element (34 1 have been proposed as gravity duals for cold atoms [24]. (See also 25 for the holographic 



description of other condensed matter systems with anisotropic scaling.) 
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This is precisely the line-element (21) of |9J. The parameters c, M in that work are related 
to the constants above by c = s/b, M = —2a 2 b. The arbitrary constant a parameterises 
the angular coordinate 0, which has periodicity 2tc. If the scaling exponent (35) is smaller 



than one, z < 1, then the line-element (37) asymptotes to AdS. If the scaling exponent (35) 
is larger than one, z > 1, then the asymptotic line-element in the limit p = £ 3 /(2r 2 ) — > oo 
takes the simple form (34) provided we choose (3 = s(£ 3 /2) z /£ 2 . If the scaling exponent (35) 
is equal to one, z — 1, then we are at the logarithmic point \fj£\ = 1 that we have neglected 
so far. 



At the logarithmic point \fj£\ = 1 we obtain two solutions. The first one is given by 



X = ( ± In p + ap + b, 0, 2/xp) 



(39) 



The angular momentum has again only a non-vanishing +- component given by J + = =F4/x. 



The result (39) leads to a spacetime that is asymptotically AdS reminiscent of the 
logarithmic mode discovered in jfjl. The other solution at the logarithmic point is given 
by 

X = (±p\np + ap + b, 0, -2 up) (40) 
The angular momentum has again only a non-vanishing +- component given by J + = 4pb, 



consistent with the result (36). The result (40) leads to a spacetime that marginally violates 



the asymptotically AdS conditions (25). It is reminiscent of the sources for operators dual 



to the logarithmic mode 26 , 27 within the logarithmic conformal field theory dual to TMG 



at the logarithmic point. Like for \p£\ ^ 1 we can fix both a and b in (39) and (40) 



by a combination of Lorentz transformations and rescalings (17) that leaves invariant the 



acceleration (30). 



4-3. Warped sector 



This sector encompasses all solutions which obey the constraints (14) and where coordinate 
vector X, velocity X and acceleration X are linearly independent. 



Imposing the constraints ( 14 ) implies 



X 



X 



XX = eij k X l X j X = €ij k X l X j X 



(41) 
(42) 



with some constant a ^ —4/£ 2 . Due to our assumption of linear independence we can write 



X = A(p) X + B(p) X + C{p) X 



(43) 



The equalities (41) imply B — C — 0. The p-derivative of a suitable combination of the 
equations of motion establishes the constraint 



XX 



XX 



(44) 
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that must be valid for all solutions. Therefore, also A must vanish in (43) and we obtain 
the condition 

X = (45) 



The general solution in this sector then is given by 

X = X(_ 2 ) p 2 + X(q) p + X(2) 



(46) 



where -X"(-2) is light-like and orthogonal to X( y, the equations of motion Q, ^ put further 



restrictions on the three constant vectors appearing in the solution (46). The Fefferman 



Graham like expansion (46) clearly is incompatible with asymptotic AdS behaviour (25). 



In fact, the solutions (46) correspond to warped (squashed or stretched) AdS geometries. 



Solutions of this type were discovered first by Nutku [28] . 

We exploit now shifts, Lorentz transformations and rescalings and obtain 



X=(p 2 - 2p, 



2^ ^ N 



(47) 



As in the previous cases the constant of motion a can be eliminated by a diffeomorphism 
together with a change of range of definition of the coordinates. The conserved angular 
momentum is given by 



3p£ 2 



(p 2 £ 2 - 9 - 8p£ 2 a + Q£ 



V, 



-2p 2 - " 



4 1)2 



9 



+ 



27 



6/i - - p 3 £ 2 - 9a + p 2 £ 2 a) (48) 
3 



The angular momentum becomes light-like for p£ = ±3 or p = 2a ± Va 2 + 9 /P. In the 
former case we obtain null warped solutions, in the latter case extremal warped black holes 
(see below). We note that the null warped solutions actually belong to the Schrodinger 
sector since X, X and X depend linearly on each other. 



To obtain the warped line-element, we take the solution for X (47) and plug it into the 
definition of the line-element (|2]), 



ds 2 = (p 2 -2p)(dx + ) 2 +(^-y 2 )(dx~) 2 -( 



App J 
p 



) dx + dx~ 



dp 2 



(49) 



4^ 2 



with M := 2p 2 - 18/ £ 2 - 4pa and J : 
x~ = a(fi leads to Nutku's solution 
A — 1 11)2 



2 pa. The coordinate transformation x + = t/a, 



28 



provided we choose a 1 



-P 2 - 

2a r 



6p 



(M 2 -9A) 



and define 



l/£ , r 2 := 2p. Other coordinate systems useful for the representation of warped AdS 
and warped AdS black holes are collected in 29 . The two zeros in the denominator of the 



dp term in the line-element (49 ) correspond to Killing horizons. An extremal Killing horizon 



emerges if the coupling constants and constant of motion are tuned as p = 2a ± ^a 2 + 9// 5 



This is precisely one of the conditions for J (48) becoming light-like that we found above. 



Warped AdS spacetimes have attracted a lot of interest recently as candidates for stable 
TMG backgrounds |29|. They arise also as string theory solutions 30 and as solutions to 



the Einstein-perfect fluid field equations 31 . 
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4-4- Generic sector 

In the Hamiltonian analysis we encountered the interesting feature that the physical 
phase space is six-dimensional in general, but only four-dimensional if the constraints 

.2 



X = XX = hold [see (14)]. These constraints are fulfilled for all sectors discussed 



so far. Note also that these sectors exhaust all solutions for which (14) holds 



The generic sector consists of all solutions violating these constraints. 

Generic: X V and/or XX ^ (50) 

It is not clear a priori if there are any solutions in the generic sector. Neither is it clear that 
this sector is empty. 

It is relatively straightforward to prove that no generic solutions exist for X that are 
finite polynomials in p. Indeed, solutions for X that are linear polynomials in p lead to 
the Einstein sector, quadratic polynomials to the warped or Schrodinger sector and certain 



higher order polynomials to the Schrodinger sector. We demonstrate in Appendix B that 
all higher order polynomial solutions belong to the Schrodinger sector. Thus, any generic 
solution must be non-polynomial in p, which makes it hard to obtain solutions by trial and 
error. In the next section we resolve this issue by providing numerical algorithms that allow 
to construct all local solutions in the generic sector. 

5. Construction of all stationary axi-symmetric solutions 

In the previous section we classified all solutions into Einstein, Schrodinger, warped and 
generic. Moreover, we constructed explicitly all local solutions in the first three sectors. 
However, we did not discuss the generic sector in detail. The purpose of this section is to fill 
this gap. Namely, we provide suitable initial data and algorithms to construct numerically 
all local solutions in the generic sector. We postpone an execution of these algorithms to 
section [6] below, where several examples will be presented. 

In the following we exploit shifts, Lorentz transformations and rescalings to simplify the 
construction of solutions as much as possible. We consider first solutions with an analytic 
centre at some point, fixed to p = 0. The solutions are parameterised by one constant of 
motion. Then we study solutions that do not necessarily have an analytic centre. We find 
that they are parameterised by three arbitrary constants of motion. In this case it is not 
possible to exploit shifts in p. Therefore, the solutions with analytic centre correspond to a 
two-dimensional subset of the three initial data. We provide simple algorithms to construct 



all local solutions and collect them in Appendix D 
5.1. Solutions with analytic centre 

We assume for the time being that there exists a centre, in the sense that X 2 = for some 
finite p. Moreover, we assume that the vector X is analytic in p at the centre. We exploit 
shifts to make X 2 = at p = 0, and Lorentz transformations to bring X into a convenient 
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form. 



X 



p=0 



(1, 0, 0) 



(51) 



Residual combinations of Lorentz transformations and rescalings can be used to simplify X 



without changing (51): we can boost in the X + X -plane and simultaneously rescale, or we 



can boost in y-direction, make a spatial rotation and then rescale. 



We study first the case where XX ^ at p = 0. In this case our Ansatz for the initial 
data consistent with the equations of motion ([7])-((9]) is 



lp=0 



X 



J 





f a 






' 3 J?t 2 


=F a \ 




±1 




j = n | 








V 2 


J 




v ° 


I 



(52) 



The evolution equations ^ do not immediately determine X, since the last term in those 
equations is multiplied by X 2 , which vanishes by our assumptions above at p = 0. However, 
taking the p-derivative of the evolution equations and using the on-shell identity 

6 



R = 2XX 



we obtain 



-ijk 



X 3 X k 



5 

2p 



3 x* 

2 



XX) Xi 



3 

" 2p 



XX) Xi 



x 2 x 



(53) 



(54) 



At p = the right hand side of (|54j) vanishes and setting the left hand side to zero determines 
X at p = by means of the matrix equation 

9 



AikX k 



"- 2 + i)* ( 



8 M (x r- 



with the matrix 



.4 



Ik 



-ijk 



X^ — XjXh 

2/i 



2/i 



at p = 



[XX)r] ik 



(55) 



(56) 



The matrix in the equations of motion (55)-(56) is invertible for XX ^ 0. Solving these 



equations establishes X and X. Taking higher p-derivatives of (54) establishes linear 
equations for higher derivatives of X at p = 0. We note that the norm of the vector 
X calculated in this way does not vanish in general at p = 0. It vanishes only if one of the 
following conditions holds 



±4 (1 



1 



) 



or 



±- 



32 



(57) 



pH 2J " ~ ~ 9 

For the same values of a it turns out that the constraint XX = holds at p = 0. For 
a = ±32/9 we obtain warped solutions, while for a = ±4 (1 — l/p 2 £ 2 ) we obtain Einstein 
solutions. 

The generic sector contains only those solutions where a is not given by the values listed 



in (57). We present in Appendix D an algorithm to construct these solutions (see Algorithm 
1). This algorithm involves the inversion of 3 x 3 matrices as the most complicated step, 
and therefore is not expensive in terms of computer time. 
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Let us now assume that XX = at p = 0. Then our Ansatz for the initial data consistent 
with the equations of motion ([7])-((9|, with no loss of generality, is 

\ 



X 



p=0 



1 ^ 




( s\ 


( 





X\ n = u 





J = n 






[a J 


V 



-a+ fa 2 - 





3 



(58) 



/ 



Here s = ±1,0 and a is a constant of motion. In the present case the matrix (56) is non 



invertible. We then take the p-derivative of equation (54) and obtain in this way an equation 
for X at p 



0: 



(A ik --XXrj ik ) X k 



l> 



9 

8/x 



9 



(X + -)X l + —XXX l -A ik X k atp = (59) 



4/x 



Solving the equations above with the Ansatz (58) turns out to be possible only if the initial 
data are further constrained. Either a = ±2/p,£ and p,i ^ ±5, in which case we find X = 
and recover solutions of the Einstein sector with light-like angular momentum; or a = 2/3, 
in which case we find XX = X =0 and recover solutions of the warped sector with 



light-like angular momentum. Thus, the Ansatz (58) does not lead to any solutions of the 
generic sector, except for a = 2/5 and \p,£\ = 5. 

Generic solutions with a = 2/5 and \/i£\ = 5 are rather special. At a certain value of p, 
which for convenience we shift to p = 0, simultaneously the following quantities vanish: 



X 2 = XX = XX = e ijk X % X 3 X 







at p = 



(60) 



However, XX ^ 0. Conversely, one can prove that all generic solutions compatible with 
these conditions require a = 2/5 and \fi£\ = 5. These solutions need the specifi cation of a 
single constant of motion, for instance the value of X~ at p = 0. We present in Appendix 
D]an algorithm to construct these solutions (see Algorithm 2). As in the previous case this 
algorithm involves the inversion of 3 x 3 matrices as the most complicated step. 



5.2. Completely generic solutions 

If there is no centre then X 2 < everywhere, and we cannot use shifts to make X light-like 
at p = 0. If the centre is not analytic we cannot use it straightforwardly as starting point 
for the time evolution. Finally, even if there is an analytic centre we may choose to start 
our evolution from a point that is not the centre. In either of these cases we have to assume 
that X is space-like at our starting point of the evolution, which we shift to p = 0. We can 
still exploit Lorentz transformations to bring X into a convenient form. 

A simple counting reveals that we have now three initial data: we need nine data to 
specify X\ p=0 , X\ p=0 and J. However, these data are constrained by (|8]) and Q. Three 
Lorentz transformations and a rescaling eliminate four additional data, so we end up with 
9 — 2 — 4 = 3 initial data. By comparison, in the case where we evolve from an analytic 
centre there are three linearly independent constraints on the initial data (|7|-d9J)- Moreover, 
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there is the constraint that X\ p= q is light-like. Again we have three Lorentz transformations 
and a rescaling. In total we have 9 — 3 — 1— 4 = 1 initial datum. An equivalent way to see 
this is as follows: there is only one residual GL(2, R) transformation that leaves invariant 
the space-like Ansatz X\ p= q = (0, 0, 1), which we can use to fix one of the components of 
X\ p=0 . So we have two initial data residing in the velocity, a, b, and one datum residing 
in the angular momentum, c + . However, there are two residual GL(2, R) transformations 
that leave invariant the light-like Ansatz X\ p=0 = (1, 0, 0), which we can use to fix two of 
the components of X\ p= q. So we have one initial datum residing in the velocity, a, and no 
datum residing in the angular momentum, as it is fully determined by ([7]). 

If the inner product between X and X is non-vanishing we can scale it to a certain 
value in this way. Then our Ansatz for the initial data, assuming b ^ 0, is 



X 



p=0 





w 



X 
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b 
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b T \ + c + 
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o =V 
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J = p 
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b 
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±1 
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— | ab 


J 



(61) 



If XX = at p = our Ansatz for the initial data, assuming b ^ 0, is 

/o\ / b\ ( b + c^ 



X 



p=0 





v 1 ; 



X 



p=0 



a 




6 ( Jp + aC+ ) 



\ 



| ab 



\ 



J 



(62) 



The initial data are parameterised by three constants a, b, c + and are consistent with the 
constraints @, Q. 

The special case 6 = requires separate discussion. We assume first a^0. If the inner 
product between X and X is non-vanishing then our Ansatz for the initial data is 



X 



p=0 





w 



X 





/ 
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= /i 
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J = n 
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±1 
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a =F | + c 
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(63) 



If XX = at p = our Ansatz for the initial data is 



\p=0 





V 1 / 



X 



p=0 



a 




/ 



V 



aPp? 

-a + c" 




/ 



(64) 



The initial data are parameterised by two constants a, c~ and are consistent with the 
constraints ([8]), (|9]). If additionally a = then there is no generic solution; instead, we 
recover solutions of the Einstein sector. 



We present in Appendix D an algorithm to construct generic solutions (see Algorithm 
3). This algorithm involves solving numerically a coupled set of three second order ordinary 
differential equations as the most complicated step. Of course, one has to check whether the 
solution obtained is really generic or belongs to some of the special cases discussed above. 
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A simple check to rule out Einstein, Schrodinger and warped solutions is to calculate X 
and XX at some point. If either of these scalar quantities is non-vanishing the solution 
necessarily is generic. If both of them vanish at one point they must vanish at all points 
and the solution is not generic. If the solution has no centre then, obviously, it cannot be 
a solution with analytic centre. If it has a centre then typically the solution will not be 
analytic there. This is so, because for some random initial data a, b, c + the subset of data 
consistent with an analytic centre is a set of measure zero. 



5.3. Summary of generic sector 

Let us now summarise the structure of the generic sector. After using rescalings and Lorentz 
transformations, the full space of solutions is parameterised essentially by three constants 
appearing in the initial data (61)-(64). A 2-dimensional surfac^] in this space of initial 



data yields the solutions with an analytic centre (52). On this surface there are curves 
corresponding to non-generic solutions, either Einstein or warped AdS. Besides the solutions 
on this surface there are two other curves in the space of initial data corresponding to 
Schrodinger solutions. For odd integer values of \ii one of the Schrodinger lines lies on the 
surface corresponding to solutions with analytic centre. To the best of our knowledge no 
solutions in the generic sector are known so farj^j In section [6] below we construct numerically 
the first examples for generic solutions. In |Appendix C we collect some interesting features 



of solutions for special values of the dimensionless coupling constant 



6. Examples 



We implement now the algorithms described in the previous section and in Appendix D 
For the generic case in subsection 5.2| our numeric algorithm essentially consists of letting 
Mathematica 33 solve the coupled system of equations ^ with initial data constrained 
by ^ and ([9]). When we solve around an analytic centre we implement equivalents of 
Algorithms 1 and 2 described in Appendix D| to obtain a Taylor series around the centre. 
Then we use this series to provide initial data with X 2 ^ close to the centre, and use the 
generic algorithm from there on. 

Then we 



First we test our algorithms and recover known solutions in section 6.1 



give examples of generic solutions with non-analytic centre in section 6.2 In section 6.3 



we 



discover soliton solutions with no centre. Finally, in section 6.4 we focus on soliton solutions 
with analytic centre. 

We plot various quantities for different solutions. Typically, we plot the components of 



Note that even if there is just one initial datum in ( 52 ), this set corresponds to a two dimensional subset of 



the initial data of (61)-(64). The reason is that one can use any value of p as starting point. Put differently, 



in (52) we fixed the shift symmetry in p which was left unfixed in (61)-(64) 



An interesting method to solve massive gravity field equations in the presence of a Killing vector was 



presented recently by Giirses 32 . It is not clear, however, how to obtain solutions of the generic sector in 
this framework. 
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X, the quadratic curvature invariant 

K = R^B!" - -R 2 = -(XX) 2 - -X 2 X 2 - -X 2 (XX) - —(X 2 ) 2 (65) 
8 2 V ; 2 4 v ; 32 v ; v ; 

and as consistency checks the Hamilton constraint G and the components of the equations 

of motion ^ 

EOM := X + — (3X x X + 2X x X) (66) 

2 2 

For solutions in the generic sector we also plot X , XX, X and the components of X. 
Some important plots are inserted in the main text, whereas others are collected in |Appendix 
[Pi after the references. 

6.1. Recovering Schrddinger and warped 

As a test of our algorithms we first recover numerically several solutions that are known 
analytically. We present one Schrddinger and one warped example. Both are obtained 
with the generic Algorithm 3. We have also recovered known solutions by making a Taylor 
expansion around an analytic centre and then starting a generic evolution from a nearby 
point, but in order to restrict the volume of the paper we do not include such an example 
in the plots. 



The plots on the left in Fig. D2 represent a numerically computed Schrddinger solution. 



For comparison, on the right the corresponding exact solution is plotted. This example fits 



in the general expression (33) with p = 0.2, £ = 1, a = 0, 6 = 5 and s = +1. We started 



the evolution from p = 1 in (33), but this point is shifted to p = in the plots. This 
means that the (non-analytic) centre is at p = — 1. As indicated by the plot of EOM, the 
numerics breaks down at this point. This happens even though the spacetime is smooth 
there. Indeed, p = — 1 does not represent a curvature singularity as is known from the 
analytic solution. The breakdown of the numerics close to the centre stems from our using 
^ as evolution equation. 



In the warped example presented in Fig. D3 similar remarks apply. However the centre 



(also here shifted to p = — 1) is analytic. This example can be obtained from (47) by setting 
p = 4, £ = 1 and a = 2. 



In all cases above the numerical results for the quadratic curvature invariant (65) match 
with high precision the exact expressions. 

K =K = ^— K -K= ~ ?2/i 2 l 2 + 243 

Schrodingcr J ^ Einstein 2^?^ warped 54^^ 

Let us now turn to our real goal in this section: describing solutions that are neither Einstein, 
nor Schrddinger nor warped. 

6.2. Generic solution with non-analytic centre — Naked singularities 

A situation that occurs frequently when one chooses random initial data is that the solution 
develops curvature singularities. These typically manifest themselves in the form of non- 



analytic centres. In Fig. D4 we show a representative example of this class of solutions. 
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For this solution we chose p — 7, £ — 1, X(Q) = (6, 5, 7), X(0) = (6, 7, 5) and J + = 8. 
With these initial data the solution develops two non-analytic centres at p w —3.3 and 
p m 1.21. That they are centres is clear from the plot of X 2 . Furthermore, the square 
invariant K diverges at these points. Thus we have a region of smooth spacetime between 
two naked curvature singularities. 

... .-2 

As is clear from the plots of XX and X this solution belongs to the generic sector. 
Not surprisingly, the numerics is not reliable close to the singular points. 



6.3. Generic solution with no centre — Solitons 

A class of interesting solutions consists of spacetimes without centre, i.e., for which X 2 < 
everywhere. Varying the initial data by trial and error, it is not hard to obtain such solutions. 



Typically, there are two asymptotic regions where the constraints (14) approach zero, and 
a central region where they deviate from zero appreciably. 

There is of course an intrinsic problem with constructing these soliton solutions 
numerically. Since we solve the equations of motion only for a limited range of p it is 
always conceivable that the solution develops pathologies for p — > ±oo. Below we remedy 
this deficiency in the most crude way: we plot a region as large as the numerics allows. This 
gives already very believable results, and allows to determine the asymptotic behaviour of 
the solutions. The next natural steps are to construct these solutions in an asymptotic 
expansion, that — as we shall see — must be non-analytic in p, and to perform a stability 
analysis. We leave these tasks for future work. 

With this disclaimer, let us now present a typical example. We choose the data p = 4, 
i = l, X(0) = (-1,7/3,-2/3), X(0) = (-1/5,-1/10,-8/3) and J+ = -11/4. In Fig. [j] 
the region near p = is plotted. X 2 has a parabolic shape and never reaches zero. The two 

... -2 

"constraints" XX and X vary substantially and take values of order one, but approach 
zero as p increases. Furthermore, K is clearly non-constant. Most of the interesting features 
in these plots are localised in a region around p ~ —0.5. 

Zooming out to \p\ ~ 10 4 produces Fig. 2 We see that asymptotically XX ~ and 

■2 

X ~ 0. We therefore expect that the solution asymptotes to a known solution in one 
of the non-generic sectors. Inspection of the components of X reveals that they approach 
constants, indicating that the solution is asymptotically warped AdS. Also K seems to 



approach the exact value K [see (67)] for warped AdS. 

While still small, already at these values of p the constraint G and equations of motion 
EOM start to deviate from zero due to numerical errors. It turns out that there is a simple 
fix for this problem. The numerical errors depend heavily on the particular Lorentz frame 
we use for our computations. Asymptotically, the preferred frame is that in which the 
acceleration X has only one non-zero component. To improve the numerical accuracy we 
therefore extract the approximate asymptotic value of X and pass to a frame where Y = 

•■2 

and X~ rj 0. Of course, we cannot put both exactly zero since X is only approximately 
null. In this new frame numerics is more stable. The whole process can then be iterated 
to get closer and closer to the optimal frame and to higher and higher values of p without 




Figure 1. The plots show a solution in the generic sector without a centre. Here [i = 4 
and 1 = 1. Evolution starts at p — where the initial values X(0) = (—1,7/3,-2/3), 
X(0) = (-1/5,-1/10,-8/3) and J+ = -11/4 were chosen. 
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Figure 2. The same solution as in Fig. [ljbut with a much larger range in p. We have let 
it evolve until the numerics becomes dubious. The facts that X approaches a constant and 
that K tends to the exact value for warped AdS at fj,£ = 4 (straight line) indicate that the 
solution is asymptotically warped AdS. 
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compromising the numerical stability. 

In this way it is possible to extend the numerics as far as necessary. We evolve until 
|p| ~ 10 12 . To understand the asymptotics of the solution we present several quantities 
plotted in Fig. [3j Three of the quantities have been multiplied with suitable functions of p 
to aid the eye in identifying the asymptotic scalings: 

(K-K),/p, XXp 3 / 2 and X 2 p 5/2 (68) 

We find that all these quantities to very high accuracy have sinusoidal behaviour in logp 
asymptotically. Therefore we are led to postulate the following asymptotic behaviour of the 
soliton solution: 

K — K ~ p~ 1//2 sin [uj log p + 81) 

XX ~ p~ 3/2 sin (w log p + 5 2 ) (69) 
X. 2 ~ p~ 5//2 sin (u log p + 5 3 ) 

for some constants u and 5{. It is clear, and also explicitly verifiable, that the subleading 
terms in the components of X also exhibit a similar functional behaviour. 



6.4- Generic solution with analytic centre — Solitons 



To exemplify the procedure of first Taylor expanding around an analytic centre and then 



using the thus obtained values as initial data we present the plots in Figs. D5 and D6 Both 
figures illustrate the same solution corresponding to a = 3.7 and the upper sign in (52), with 
p = 5 and 1=1. The same qualitative behaviour is obtained for the special case obtained 
when setting a = 2/5 in ([58|)f| 



Fig. |D5| shows the region close to the centre at p = 0. Here we include both the Taylor 

expansion itself and the solution obtained from using the values at p = —0.1 in the Taylor 

expansion as initial values. It is clear that the solution is in the generic sector since neither 
•■ 2 

XX nor X vanishes. 



In Fig. D6 the range of p is extended to negative values until the numerics starts to be 

... .-2 

questionable. Again, the two quantities XX and X rapidly approach zero as p — > —00, 
and X approaches constants. 

Also for this solution it is possible to deduce the asymptotic behaviour by going to a 
different Lorentz frame. The result is very similar to the example in subsection 6.3| with 



exactly the same scalings (69) (but of course with different u and 5i). Note that this implies 



that the half-integer exponents of p in the asymptotic expansion (69) are not functions of pi 



This can also be verified by choosing non-integer values for pi. The solution with analytic 
centre thus is a soliton that asymptotes to warped AdS. We give the corresponding plots in 
Fig-i 



5 In this special case we use first Algorithm 2 described in Appendix D for the Taylor expansion. Then we 
extract initial values at some value of p that we feed into Algorithm 3, exactly as described in the text. 
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Figure 3. This figure shows the same solution as in Fig. [T]but Lorentz rotated to stabilise 
the numerics. We see that the numerical errors are small as compared to Fig. [2] Note 
that the scale on the p-axis is logarithmic and that we have multiplied three quantities by 
half-integer powers of p. Asymptotically these quantities exhibit sinusoidal dependence on 
log p. 



In this example we choose to include a plot (top-right) that reveals the subleading 
behaviour of X + : 

X + = Ap 2 + Bp 3/2 sin (wlogp + cf) + ... (70) 
Note in particular that the next-to-leading contribution grows faster than p which makes 



the solution incompatible with Compere-Detournay boundary conditions 34 
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Figure 4. The same solution as in Fig. |D6| is plotted in a more favourable Lorentz frame. 
We observe the same oscillating behaviour as in the case without a centre. For this case 
we also include a plot of \p\ 3 ^ 2 d p (X + / p 2 ) from which the sub-leading behaviour of X + 
can be deduced. Note that the errors are small and that we plot against \p\ = —p, so the 
asymptotic region is to the right in the plots. 



6.5. Discussion and open questions 



The soliton solutions discovered in section 6^3, OA asymptote to warped AdS in the sense 
explained above. At first glance this appears to contradict the no-go result of [llj, whose 
analysis in Appendix E can be generalised straightforwardly to non-null- warped AdS. 
However, this no-go result relies on the assumption that the solution can be represented as 
an asymptotic expansion in negative integer powers of p. Our soliton solutions are capable 
to circumvent the no-go result since the subleading terms in our asymptotic expressions are 



not analytic in the coordinate p (69), (70). 



The existence of these soliton solutions raises several further questions that are beyond 



All stationary axi- symmetric local solutions of topologically massive gravity 



24 



the scope of the present paper. Let us end this section by listing a few of them. 



Are there consistent boundary conditions that can accommodate the asymptotics (69) 
(70)? If so what is the corresponding asymptotic symmetry group? 



Which boundary terms are needed to yield a vanishing first variation of the action? 

Are the solutions stable with respect to perturbations that do not respect axi-symmetry 
and/or stationarity? 

Are there any exact solutions of this kind? As the asymptotic functional form is 



established (69), (70) it is conceivable that exact solutions can be constructed. 
Can these solutions be interpreted as finite energy excitations around warped AdS? If 
this question can be answered affirmative then the soliton solutions found above could 
have natural interpretations as geons or warped black holes with topologically massive 
graviton hair. 



• What is the precise range of initial data (61)-(64) for which we obtain smooth soliton 
solutions as opposed to naked singularities? 

• Are there soliton solutions that asymptote to AdS or Schrodinger? 

• Are there kink solutions that interpolate between two different non-generic 
backgrounds, say between AdS and warped AdS? The existence of such solutions could 
be used to address the stability issues]^] 

This concludes our presentation of selected numerical solutions of TMM. In summary 
we have found new soliton solutions in TMG, unravelled their asymptotic warped AdS 
behaviour and explained why the existence of these solutions does not contradict previous 
no-go results. We have also raised several questions to be addressed in future work. In 
the next section we generalise our discussion to non-negative values of the cosmological 
constant. 



7. Generalisations 

The discussion so far can easily be generalised to vanishing {t — > oo) or positive (£ 2 < 0) 
cosmological constant. We start with the latter case. 

We classify again the solutions into Einstein (X = 0), "Schrodinger" (X ^ 0, 
X = a(p)X + P(p)X), warped (X 2 = XX = 0, e ijk X i X^X k ^ 0) and generic (X 2 ^ 
and/or XX ^ 0). 



7.1. Positive cosmological constant 

The essential difference to the AdS case is the sign change of the constant term in the 
Hamilton constraint Q. 

6 Kink solutions were found in a related 3-dimensional theory, new massive gravity, by Oliva, Tempo and 
Troncoso [49 . Evidence for kink solutions in TMG was reported by Lashkari and Maloney studying the 
Ricci-Cotton flow in TMM [50]. 
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Einstein The Einstein sector persists, but there are no longer black hole solutions. All 
solutions are locally dS. 

Schrddinger The Schrddinger sector does not persist. In other words, there are no solutions 
where the acceleration X depends linearly on velocity X and coordinate vector X. There 
is a simple geometric reason for this. Namely, the Hamilton constraint ^ requires X to be 



time-like. On the other hand, the constraints (14) must hold and they imply that X is null 



or space-like. Thus, there are no Schrodinger solutions for positive cosmological constant. 



Warped The warped sector persists. In fact, one can take the warped AdS solution (47) 
and just insert negative values of £ 2 to obtain warped dS solutions. 

Generic One can use here algorithms analog to the ones described in section [5] and 



Appendix D 



7.2. Vanishing cosmological constant 

The essential difference to the AdS case is that the constant term vanishes in the Hamilton 
constraint Q. 

Einstein Also here the Einstein sector persists, but no black hole solutions arise. All 
solutions are locally Minkowski. 

Schrddinger The "Schrddinger sector" is modified appreciably as the Hamilton constraint 
Q implies that X is null. The constraints (14) again hold automatically. Up to Lorentz 



transformations, rescalings and shifts the unique solution is given by 

X = (se-» p + ap + b, 0, ±l) (71) 
which depends on two integration constants a, b and a sign s — ±. Analog to the AdS- 



Schrodinger case (33) we can still exploit GL(2, R) transformations to fix a,b conveniently. 
The solution (71) is well-known [35], and it does not exhibit asymptotic Schrodinger 
behaviour. The angular momentum has only one non-vanishing component J + = a. 

Warped The warped sector leads to Vuorio-type of solutions [36] . The solution is obtained 



from (47) in the limit £ — > oo. 



Generic One can use here algorithms analog to the ones described in section [5] and 



Appendix D In fact, a numerical study of soliton solutions for this case was performed 



in 



35 with results similar to those obtained in this work for the cosmological case: Some 



initial data were found to lead to singularities, while other data lead to smooth solutions. 
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7.3. Generalisations beyond topologically massive gravity 

Clement's Ansatz for the line-element ^ can be used also for other 3-dimensional theories 



of gravity, like new massive gravity (NMG) [37] generalised massive gravity 39 
supergravity 



40 , higher order massive gravity |41 , extended new massive gravity 42 



massive 
or 



Born-Infeld massive gravity 43 . This could allow to construct stationary axi-symmetric 
solutions for these novel theories along the lines of our paper: first reduce the gravitational 
action to a particle mechanics action (see section [2]), next perform a Hamiltonian analysis 
(see section [3] and Appendix A ) , then discuss various simple sectors of the theory — analog 



to the Einstein, Schrodingeijj and warped sectors discussed in the present work (see section 
[1]) — and finally provide algorithms to construct generic solutions numerically (see sections 
[5j [6] and Appendix D). It would be interesting to construct soliton solutions for these 3- 
dimensional theories, to study their stability properties and to address the questions raised 
at the end of section 16.51 
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Appendix A. Poisson brackets and Hamilton equations of motion 



The non- vanishing Poisson brackets between all twelve constraints (12) read as follows. 

1 



{ P *,G} = -n^o 



{vlx^j} = -Xi - — e ijk x j y k ( A - la ) 



{G,x j t> j } = -y i t> i 
3 



{G l n i } = $ j »0 {G,$ j } = t> j 

1 13 
{Ui, Uj} = — e ijk x k {IT, $, } = r} i:j + — e ijk y k {Hi, x 3 ^j} = — e ijk :r : 

2 

$ 3 } = e ijk z k tfVj} = -% 

I 1 



Here are some additional Poisson brackets of interest. 
{plV^j} = -Vi 



{plz^ ] } = -2z l 
{n h y j ^j} = - Zi 



(A.lb) 
(A.lc) 

(A.ld) 



(A.le) 
(A.lf) 



7 Indeed, Clement has used this Ansatz recently [38| and found warped black hole solutions in NMG. 

8 For NMG, this sector contains both Schrodinger and Lifshitz solutions. 44 45 . 
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All other Poisson brackets between the constraints (12) and the quantities y l ^i, z l ^i vanish 



The extended Hamiltonian ( 13 ) with arbitrary A* commutes weakly with all constraints, 
except with the following ones: 

{i? ext ,$i} « e ^i--e m Xx k ) {# cxt ,x^ 4 } « -e (y% + X' Xi + ^- e ijk A Vj/*) (A.lg) 

Setting to zero the right hand sides in these equations establishes results for the Lagrange 
multipliers A*. With the Ansatz 



A 1 = X x x % + X y y l + X z z l 



(A.2a) 



we obtain 

.,2 



A* — Cij-fc xyz fc 



y% (A.2b) 



1 

Note that by assumption of regularity x 2 ^ 0, so X x is well-defined. If x, y, z are linearly 
independent then also X y and X z are well-defined. If in addition the conditions (14) are 
fulfilled then A* = and the extended Hamiltonian H ext reduces to the canonical one 



H = eG. The Poisson brackets between the canonical Hamiltonian and the conditions (14) 



either vanish or yield something proportional to these conditions, see (A. If). Therefore, if 



the conditions (14) both hold initially they must hold for all times. If x,y,z are linearly 



dependent then the constraints (|14|) are fulfilled automatically, as discussed in the main text 
in section I 



The Hamilton equations of motion generated by the extended Hamiltonian ( 13 ) are 

(A.3a) 



x 







e z 



eX l 







Vi 

Pi 



3 ~M 



£iik V Z 

2fJL 3 



e 

'2/1 



^■ijk % Z 







(A.3b) 
(A.3c) 
(A.3d) 



We see that the choice (13) for the extended Hamiltonian automatically leads to the desired 



gauge condition e = const. With no loss of generality we set e = 1 throughout in this work. 



Appendix B. Proof that all higher order polynomial solutions are Schrodinger 

To show the absence of generic solutions that are polynomial in the radial coordinate p we 
insert the Taylor series X = 52 n=0 X ( n ) p n into the equations of motion (|5| and assume 
N > 3, since for iV < 2 we recover non-generic solutions. 



N N N 

<n - l)X (n)l p- 2 = -— e ijk (s £ £ nm ^ m ~ l ) X U X U ^ 

n=2 n=l m=2 



71=1 771=2 

N N 



(B.l) 



n=0 m=3 
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We then solve (B.l) order by order in p, starting with the highest power p 2N 3 . We obtain 



the following chain of conditions: 

2 v " 3 : holds trivially 



P 

2JV-4 . -yj yk _ n 



P 



2N-5 



€iik X 



ijk^{N-2)X(N) 







n 2Af-6 . e ( yj Y k -L 2 ^ 4 Y k \ — n 

P ■ € ijk{^( N - 3 )^(N)^ 7tfj~ ^(N-2)^(N-1)) - {> 



(B.2) 
(B.3) 
(B.4) 

(B.5) 



and similar relations down to the power of p N 1 . These relations involve only the right 



hand side of the equations of motion (B.l) and taken together they imply a chain of 
proportionalities 

X (n) oc X {N) Vn > 2 (B.6) 



Apart from the proportionality constants in (B.6) only three coefficient vectors are free 



at this stage: -X"(o), -^(l) an d -X"(JV)- We then continue to solve (B.l) order by order in 
p, starting with the power p N ~ 2 down to the power p°. We obtain the following chain of 
conditions: 

2N- 1 



P 



N-2 



JV-3 



(N)i 



2/i 



c YJ Y K 



(B.7) 



*{N-i)i = ~ e ijk ((2iV - 3)X( 1) X^ N _ 1) + 2NXl Q) X* N) ) (B.8) 



P 



o . 



X { 



(2)i ~ ~2p tijk V X fa X W + 6X ?o) X (3) y 
These relations taken together imply a chain of inner products 

X(i)X( n ) = X( n )X( n ) = Vn > 2 
and the vanishing of all triple products 



e ijk ^(O)^l)^(n) 







Vn 



(B.9) 



(B.10) 



(B.ll) 



Therefore, the acceleration X must depend linearly on the vector X and its velocity X. 
This is the defining feature of a Schrodinger solution, so all polynomial solutions of at least 
third order in p are Schrodinger solutions. This is what we wanted to prove. 



Appendix C. Solutions for special values of the coupling constant 

In this appendix we collect some interesting features of solutions for special values of the 
dimensionless coupling constant \p£\. We always keep I finite and positive, even when taking 
certain singular limits. We start with the (singular) limit \pi\ — » 0. Then we address the 
logarithmic point = 1), the null warped point = 3) and the special point = 5. 
Finally we discuss the (singular) limit \p£\ — > oo. 
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Appendix C.l. Chern-Simons point 







In this (singular) limit only the gravitational Chern-Simons term remains in the action 
([I]). The limiting theory does not exhibit any propagating physical degree of freedom. The 
equations of motion require the Cotton tensor to vanish. Thus any solution of the limiting 
theory that is also a solution of TMG must belong to the Einstein sector. We shall exhibit 
this feature explicitly below. 

All solutions of the limiting theory are known in closed form. They can be obtained as 
follows. Take the gravitational Chern-Simons term and perform a Kaluza-Klein reduction 



to two dimensions 46 along the angular Killing direction. 



— CL a CLf3 
-dp 



1 



(C.l) 



Here g a p is a 2-dimensional metric and a a a 2-dimensional gauge field. The scalar field that 
appears in the Kaluza-Klein Ansatz (C.l) does not play any role in the reduced theory 
and is set to unity for the time being, tp — 1. This reduction yields a specific Maxwell- 
dilaton gravity model, the solutions of which can be obtained locally and globally in closed 
form [47]. There are some isolated solutions, so-called constant dilaton vacua, and a set of 
generic solutions. 



In an Eddington-Finkelstein patch the generic solution reads 47 

1 



dsi 



9a/3 



dx a dx 13 = 2 du dp + (2C 



Ip 2 + \p 4 )^ 2 



P °a 



(C.2) 



The quantities c,C are constants of motion and have a 2-dimensional interpretation as 
conserved [/(l)-charge and conserved mass, respectively. For BPS saturation 48 , c 2 = 8C, 
the solution (C.2) is the kink solution discovered in [46). Plugging the 2-dimensional solution 
(C.2) into the Kaluza-Klein Ansatz ( |C.1[ ) yields the 3-dimensional line-element 



ds 2 



2 du dp + (2C - - p 2 ) du 2 - p 2 du d<p - d<p 2 



(C.3) 



The geometry defined by the line-element (C.3) has vanishing Cotton-tensor and thus is 



conformally flat, but it has non-vanishing Ricci-tensor. The corresponding Ricci-scalar is 
not constant, R = c 



2 p 2 . Thus the solution (C.3) cannot arise directly as the smooth 



limit of some TMG solutions for finite values of \p£\. However, for non- vanishing mass, 
C 7^ 0, we can multiply the solution (C.3) with the conformal factor (p = 2C£ 2 /p 2 so that 
the transformed Ricci-scalar is constant and negative, R = —6/£ 2 . This solution then arises 
as the smooth limit of Einstein sector TMG solutions for finite values of \p,£\. 

The 2-dimensional equations of motion exhibit a Z2-symmetry a a — > —a a [46]. One of 
the isolated solutions preserves this Z 2 -symmetry 

c 



ds 2 = 2 du dp + (A + Bp - - p 2 ) du 2 - d0 2 



(C.4) 



while the others break it 



ds 2 = 2 du dp + (A + Bp) du 2 ± 2^~cp du d(j) - d0 2 



(C.5) 
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In all cases A, B,c are constants of motion, the Ricci-scalar is constant and the Cotton- 
tensor vanishes. For c = A/£ 2 the Z 2 -symmetry breaking solutions (C.5) directly arise 



as the smooth limit of Einstein sector TMG solutions for finite values of \fi£\. For other 



values of c 7^ we have to multiply the line-element (C.5) by a constant conformal factor 
if = c£ 2 /4 in order to obtain the appropriate transformed Ricci-scalar R = —6/£ 2 . Similarly, 



we have to multiply the Z^-symmetry preserving line-element (C.4) by a conformal factor 
if = £ 2 (A 2 + 2Bc)/(2c(p — A/c) 2 ) in order to obtain the appropriate transformed Ricci-scalar 
R = —6/£ 2 . This is possible as long as c ^ and A 2 + 2Bc ^ 0. 

In conclusion, the solutions of TMG at fj£ = form classes of conformally related 
metrics. In nearly all such classes there is a representative that can arise as the smooth 
limit of Einstein sector TMG solutions for finite values of \/j,£\. 



Appendix C.2. Logarithmic point 



At this point Einstein solutions necessarily have light-like angular momentum (21). 



Schrodinger solutions have a logarithmic term in p and can either be asymptotically AdS 



(39) [compare with (25)] or marginally violate the condition of asymptotically AdS (40). 
Neither warped nor generic solutions do exhibit special features at this point. 



Appendix C.3. Null warped point = 3 

Schrodinger solutions that are not asymptotically AdS have a scaling exponent z 



2 (35). 



The coordinate vector X is quadratic in the coordinate p (33) for these solutions. The 



angular momentum is light-like(36) [see also (48)]. The null warped solution coincides with 
the Schrodinger solution. This is the most distinguishing feature of the null warped point. 
Neither Einstein nor generic solutions do exhibit special features at this point. 



Appendix C.4- Special point 

Generic solutions with analytic centre can have peculiar properties. Namely, the conditions 
(60) simultaneously can hold at p = 0. We stress that this is the only value for the 



dimensionless coupling constant \p,£\ where generic solutions of this type can exist. For 



any other value of \fi£\ the constraints (60) imply the conditions (14) and thus do not lead 



to solutions of the generic sector. The algorithm to construct these solutions differs from 
the algorithms that are used at any other value of the coupling constant, see Algorithm 2 
in |Appendix D Neither Einstein, nor Schrodinger nor warped solutions do exhibit special 
features at this point. 



Appendix C.5. Einstein point \fj£\ — > oo 

In this (singular) limit only the Einstein-Hilbert term and the cosmological constant remain 
in the action Q. The limiting theory does not exhibit any propagating physical degree of 
freedom. Clearly, all solutions in this limit are in the Einstein sector and therefore locally 
and asymptotically AdS. 
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Appendix D. Algorithms 



In Fig. Dl we depict a flow chart for the decision which of the algorithms below should be 
used (if any), given some initial data X, X at p = and the conserved angular momentum 



J. For instance, if the constraints (14) both hold then no algorithm is needed since we have 



a non-generic solution, that is, either Einstein or Schrodinger or warped. These solutions 



are all known analytically. If the constraints ( 14 ) do not hold then for random initial data 



one can use the Algorithm 3 below. If the initial data are fine-tuned so that a regular 
centre emerges and the starting point of the evolution coincides with the centre then one 
has to use Algorithms 1 or 2, depending on the outcome of the flow chart. The matching 
procedure indicated by the dashed arrows works as follows: we follow Algorithms 1 or 2 
up to the required level of iteration, say, hundred steps. In this way we know the value of 
the first hundred derivatives of X at p = and thus we can construct the solution at some 
finite value of p for a given accuracy, say, at p = 0.01. We then take the values of X, X at 
p = 0.01 and the values of J as new initial data and feed them into Algorithm 3. 

In the algorithms below we assume that Lorentz transformations and rescalings have 
been exploited as explained in the main text in order to reduce the amount of independent 



initial data. However, we stress that the flow chart in Fig. |D1| works for any set of initial 
data, and that the algorithms below easily can be adapted to initial data that has not been 
reduced by these transformations. 



Algorithm 1 — X 2 = 0, XX ^ at p = 



Step 1: Initial data One has to choose the upper or lower sign in the Ansatz (52) and 



a value for a that differs from the values listed in (57) and is compatible with X being 
spacelike for p > 0. 



±/i < 



^ ±4 (l-7^) 



a ^ ± 
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p, 2 P J ~ ' ~ 9 

The initial data (52) by construction are compatible with the constraints ([7])- ([9]). 



(D.l) 



Step 2: Solving for X Since the matrix (56) is invertible we obtain a unique solution for 



X at p = by solving the linear equations (55). 



Step 3: Solving for X Similarly, we obtain a unique solution for X at p = by solving 



the linear equations (59). Note that the matrix appearing on the left hand side of (59) 



necessarily is invertible. 



Step n: Solving for d n X / dp n By taking the n-th p-derivative of (54) and evaluating it at 



(n) 

p = we obtain three linear equations that allow to determine X at p = 0. The crucial 
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Figure Dl. Flow chart for the decision which algorithm to use starts in upper left 
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input here is that all the matrices appearing in this construction are invertible, 



A 



(n) 
ik 



A ik - ^ — ^ XXr] ik det A {n) ± V n > 2 



(D.2) 



with Ajfc as defined in (56|). We obtain the evolution equation 
9 d n ~ 2 



, ( n ) Q H™- 2 / n A \ n ~ 2 Am-1 , (n-m) 

A ( ^' = | i( ^((^ + ^)-E < ^l(^-" > »^) „>8 (D.3) 



The right hand side of (D.3) contains X and its first n — 1 derivatives, all of which can be 



determined iteratively by using (D.3) repeatedly. The matrix on the left hand side of (D.3) 



is sensitive to the initial data (52), only. 



. .. (n) 

Final Step: Summing the Taylor series Given the solutions for X, X } X, ... X at p = 
for some finite n the solution for X is approximately given by 



X^X\ n + pX\ n 

I p=0 r I p=0 



p x\ 



+ ^T X n 



(D.4) 



provided that p is sufficiently small. At some finite value of p one can then match to the 
algorithm for solutions with no analytic centre described below to evolve to arbitrary values 
of p. 



* * * 



Algorithm 2 — X 2 = XX = at p = 



Step 1: Initial data This algorithm only is applicable when \p£\ = 5. One has to choose 



s = ±1, in the Ansatz (58) and a = 2/5. A convenient choice for the free initial datum is 
the value of X~ at p = 0. 



Step 2: Consistency of X The condition (55) simplifies to A^X^ 1 = 0, where the matrix 



Aik now has rank 1 only. This condition then establishes X =0. 



Step 3: Solving for X Taking the p-derivative of (54) yields 



with the invertible matrix 

Bik = 



BikX k — AikX k 



6 19 • • 5 

Vik + ~r~ XiXk + — XiXk — eijk X 3 
pi 1 4p 2p 



(D.5) 



(D.6) 



The quadratic set of equations (D.5) can be solved uniquely for X k in terms of the initial 
datum X~ . It is better to solve first the ^-component of (D.6) for X Y and then the minus 



component of (D.6) for X + . In this way only linear equations have to be solved. The plus 



component of (D.6) is an identity so that the solution for X necessarily is compatible with 



the consistency requirement of Step 2 above. 
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Step n + 1: Solving for d n X / dp n By taking the n-th p-derivative of (54) and evaluating it 



(«) 

at p = we obtain three linear equations that allow to determine X at p = 0. The crucial 
input here is that all the matrices appearing in this construction are invertible, 



B 



(ft) 
if. 



9 • • 5 

~r~ XiX k + — XiX k 
4// 2p 



E4 m) 

m=2 



det 5 (n) ^ V n > 3 



(D.7) 



where the matrices are defined in (D.2) below. We obtain the evolution equation 



B^X K + A lk X 



(ft) 



(n+l) 



9 dn_1 f(x 2 + 4 V* 



n-2 



9 (ft) . . 
— XXXi 
4p 

n-1 



5 



2p 



XXX, 



-E 



; dp" 



(n-m) ft-1 (n) 

A ( r m) ^ fc ) + E A ( r^ fc " > 3 (D.8) 

m=2 



The right hand side of the evolution equation (D.8) contains up to n — 1 derivatives of X 



(note in particular that all n-derivative terms cancel). The left hand side of the evolution 

(ft) 



equation (D.8) does contain not only X, but the minus component of this equation also is 



(ft+i) 



sensitive to X . This component can then be determined as follows: we solve (D.8) for 



(n) 

X in terms of lower derivatives (which are known at this Step) and the minus component 

(ft+i) 

of X . Since the minus component of X is determined from Step n already, we obtain in 

(ft+i) 

this way one equation for the unknown minus component of X . 



Final Step: Summing the Taylor series This step is exactly the same as in Algorithm 1 
above. 



* * * 



Algorithm 3 — X 2 ^0atp = 

Step 1: Initial data One chooses three constants a,b,c + and starts with the initial data 
(61) or (62) if b ^ 0. By construction they are compatible with the constraints (Js]) - (J9j) - If 
b = but a ^ then one takes instead the initial data (63) or (64). If a = b = then one 



does not obtain a solution of the generic sector but rather an Einstein solution. 



Final step: Time evolution With these initial data one solves numerically the coupled 
set of second order ordinary differential equations (J7|). This can be done straightforwardly 
with some computer algebra system. As long as X 2 ^ the configurations obtained in 
this way will remain regular. If there is a value for p where X 2 — > then numerics will 
break down eventually and one cannot follow the evolution beyond this point. But such a 
locus anyhow is a boundary of our coordinate patch. It is worthwhile emphasising that any 
solution obtained in this way automatically solves the original equations of motion ^ and 
the Hamilton constraint 
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Schrodinger numerically 
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Figure D2. Plotted side by side are X, K, G and EOM for a Schrodinger solution. The 
left column of plots was obtained numerically, whereas the right column corresponds to the 
exact solution. This particular solution corresponds to the values p = 0.2, I = 1, a = 0, 



6 = 5 and s = +1 in (33 1. We start evolving from p = 1, which is shifted to p = in the 
plots. 
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Warped AdS numerically 
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Figure D3. Plotted side by side are X, K, G and EOM for a warped AdS solution. The 
left column of plots was obtained numerically, whereas the right column corresponds to the 
exact solution. This particular solution corresponds to the values p, = 4, £ = 1 and a = 2 



in (47). We start evolving from p — 1, which is shifted to p = in the plots. 
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Naked singularity solution 




Figure D4. The plots show a numerical solution in the generic sector. It has two non- 
analytic centres at p « —3.3 and p ~ 1.21 respectively. As can be deduced from the plot 
of K both centres correspond to curvature singularities. Note that the numerics breaks 
down close to the singularities. The solution corresponds to initial values X(Q) = (6, 5, 7), 
X(0) = (6, 7, 5) and J+ = 8. We choose fj, = 7 and 1=1. 
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Figure D5. The figure shows a solution in the generic sector with an analytic centre at 
p — 0. In the plots both the Taylor expansion and the corresponding extension by algorithm 
3 are shown. We have chosen fi = 5 and £ = 1, and the solution corresponds to thee upper 



sign in (52 1 and a — 3.7 which is in between the Einstein and warped values. 
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Figure D6. The same solution as in Fig. D5 but with a wider range of p. Asymptotically 
the solution looks like warped AdS. The straight line in the plot of K is the exact value for 
warped AdS. It is clear from the plots of G and EOM that the numerics are problematic 
for large negative p. 



